Light Curves Research Articles

Long-term stellar activity of M dwarfs. A combined K2 and TESS study of two early M-type stars

Studies of the rotation and activity of M type stars are essential in order to enhance our understanding of stellar dynamos and angular momentum evolution. Using the outstanding photometric capabilities of space telescopes, even those rotation signals with low amplitudes can now be investigated in previously unrivaled detail. By combining data of K2 and the TESS prime mission, the star spot activity of M dwarfs can be monitored over a timescale of half a decade. In the framework of our study on the rotation--activity relation for bright and nearby M dwarfs, we also aim to investigate their long-term activity. While K2 was observing fields distributed around the ecliptic plane, the TESS prime mission was oriented along a line of ecliptic longitude, with one camera centered on an ecliptic pole. Due to these different observing strategies, the overlap between K2 and the TESS prime mission is marginal. However, 45 stars from our sample were observed with both missions, and of these, two early M-type stars that fulfill our selection criteria, EPIC\,202059229 and EPIC\,245919787, were analyzed in more detail. We find that, for both stars, the rotation period did not change, while the rotational phase did change for EPIC\,245919787 by sim 0.2. The amplitude of the spot-induced variability changed for both stars but more significantly for EPIC\,245919787. By comparing the cumulative flare frequency distributions, we find that the flare activity for EPIC\,202059229 is unchanged, while it slightly changes for EPIC\,245919787 between the K2 and TESS epochs. Using a combination of light curves from K2 and TESS that span a baseline of up to 4.5 years, we are able to measure significant differential rotation for EPIC\,245919787. Furthermore, we show that combining missions like K2 and TESS is a promising method for detecting stellar activity cycles.

Exoplanet Transit Spectroscopy with JWST NIRSpec: Diagnostics and Homogeneous Case Study of WASP-39 b

Abstract The JWST has ushered in a new era of exoplanet transit spectroscopy. Among the JWST instruments, the Near-Infrared Spectrograph (NIRSpec) has the most extensive set of configurations for exoplanet time series observations. The NIRSpec Prism and G395H grating represent two extremes in NIRSpec instrument modes, with the Prism spanning a wider spectral range (0.6-5.3 μ m) at lower resolution (R∼100) compared to G395H (2.87-5.14 μ m; R∼2700). In this work, we develop a new data reduction framework, JexoPipe, to conduct a homogeneous assessment of the two NIRSpec modes for exoplanet spectroscopy. We use observations of the hot Saturn WASP-39 b obtained as part of the JWST Transiting Exoplanets ERS program to assess the spectral quality and stability between the two instrument modes at different epochs. We explore the noise sources, effect of saturation, and offsets in transmission spectra between the different instrument modes and also between the two G395H NRS detectors. We find an inter-detector offset in G395H of ∼ 40-50 ppm, consistent with recent studies. We find evidence for correlated noise in the Prism white light curve. We find the G395H spectrum to be of higher precision compared to the Prism at the same resolution. We also compare the JexoPipe spectra with those reported from other pipelines. Our work underscores the need for robust assessment of instrument performance and identification of optimal practices for JWST data reduction and analyses.

The structure and evolution of relativistic jetted blast waves

Abstract We study, analytically and numerically, the structure and evolution of relativistic jetted blast waves that propagate in uniform media, such as those that generate afterglows of gamma-ray bursts. Similar to previous studies, we find that the evolution can be divided into two parts: (i) a pre-spreading phase, in which the jet core angle is roughly constant, θc, 0, and the shock Lorentz factor along the axis, Γa, evolves as a part of the Blandford-Mckee solution, and (ii) a spreading phase, in which Γa drops exponentially with the radius and the core angle, θc, grows rapidly. Nevertheless, the jet remains collimated during the relativistic phase, where $\theta _c(\Gamma _a\beta _a=1)\simeq 0.4\theta _{c,0}^{1/3}$. The transition between the phases occurs when $\Gamma _a\simeq 0.2\theta _{c,0}^{-1}$. We find that the ”wings” of jets with initial ”narrow” structure ($\frac{d \log \, E_{iso}}{d\log \, \theta }<-3$ outside of the core, where Eiso is isotropic equivalent energy), start evolving during the pre-spreading phase. By the spreading phase these jets evolve to a self-similar profile, which is independent of the initial structure, where in the wings Γ(θ)∝θ−1.5 and Eiso(θ)∝θ−2.6. Jets with initial ”wide” structure roughly keep their initial profile during their entire evolution. We provide analytic description of the jet lateral profile evolution for a range of initial structures, as well as the evolution of Γa and θc. For off-axis GRBs, we present a relation between the initial jet structure and the light curve rising phase. Applying our model to GW170817, we find that initially the jet had $\theta _{c,0}=0.4-4.5~\deg$ and wings consistent with Eiso∝θ−3 − θ−4.

An analytic description of electron thermalization in kilonovae ejecta

Abstract A simple analytic description is provided of the rate of energy deposition by β-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters- initial entropy, electron fraction {s0, Ye} and scaled (time-independent) density ρt3. The formulae are derived using detailed numerical calculations following the time-dependent composition and β-decay emission spectra (including the effect of delayed deposition). The deposition efficiency depends mainly on ρt3 and only weakly on {s0, Ye}. The time te at which the ratio between the rates of electron energy deposition and energy production drops to 1 − e−1, is given by $t_e=t_{0e}\Big (\frac{\rho t^3}{0.5(\rho t^3)_0}\Big )^a$, where $(\rho t^3)_0=\frac{0.05\, {\rm M}_{\odot }}{4\pi (0.2c)^3}$, t0e(s0, Ye) ≈ 17 days and 0.4 ≤ a(s0, Ye) ≤ 0.5. The fractional uncertainty in te due to nuclear physics uncertainties is $\approx 10~{{\%}}$. The result a ≤ 0.5 reflects the fact that the characteristic β-decay electron energies do not decrease with time (largely due to ‘inverted decay chains’ in which a slowly-decaying isotope decays to a rapidly-decaying isotope with higher end-point energy). We provide an analytic approximation for the time-dependent electron energy deposition rate, reproducing the numerical results to better than 50% (typically $<30~{{\%}}$, well within the energy production rate uncertainty due to nuclear physics uncertainties) over a 3-4 orders-of-magnitude deposition rate decrease with time. Our results may be easily incorporated in calculations of kilonovae light curves (with general density and composition structures), eliminating the need to numerically follow the time-dependent electron spectra. Identifying te, e.g. in the bolometric light curve, will constrain the ejecta density distribution.

H.E.S.S. observations of the 2021 periastron passage of PSR B1259-63/LS 2883

is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star (LS 2883). At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission, for which the temporal and spectral properties of this emission are, for now, poorly understood. In this regard, very high-energy (VHE) emission is especially useful to study and constrain radiation processes and particle acceleration in the system. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of approximately 100 hours of data taken over five months, from p -24$ days to $t_ p +127$ days around the system's 2021 periastron passage (where $t_ p $ is the time of periastron). We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent overall with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve -- although sparsely sampled in this time period -- we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index $ stat sys $, a value consistent with the spectral index obtained from the analysis of data collected with H.E.S.S. during the previous observations of the source. We report spectral variability with a difference of $ stat sys $ at $95 <!PCT!>$ confidence intervals, between sub-periods of the 2021 dataset. We also detail our investigation into X-ray/TeV and GeV/TeV flux correlations in the 2021 periastron passage. We find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after $t_ p +25$ days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from $ t_p-23$ days to $ t_p+126$ days. This suggests that the GeV and TeV emission originate from different electron populations.

A study of two young multipolar planetary nebulae: Hen 2-73 and Hen 2-96

We perform an infrared (IR) spectral and visible morphological study of two young planetary nebulae (YPNe) Hen 2-73 and Hen 2-96 using archival Spitzer Space Telescope and Hubble Space Telescope (HST) observations to understand their dust properties and nebular structures. High-resolution HST images of these nebulae show several bipolar lobes and ionised tori in the central regions of both objects. The presence of these multi-lobe structures suggests that the formation process of these nebulae is complex. To search for a possible link between the central sources and multipolar appearances of these objects, the Transiting Exoplanet Survey Satellite (TESS) observations are used to examine whether their central stars (CSs) exhibit periodic photometric variations. In the TESS observations, the CS light curve of Hen 2-96 shows a photometric variation with a period of 2.23 hr. The IR spectra of these two YPNe suggest that the nebulae have mixed dust environments, which are associated with the presence of dense tori created by central binary interactions in these objects. Two three-dimensional models are constructed to study the complex nebular structures of the YPNe. These simulations suggest that the number of multipolar YPNe may be larger than observed. In addition, we analyse the spectral energy distributions of these nebulae to study their gas, dust, and photospheric components.

Unveiling the periodic variability patterns of the X-ray emission from the blazar PG 1553+113

The search for periodicity in the multiwavelength, highly variable emission of blazars is a key feature to understanding dynamical processes at work in this class of active galactic nuclei. The blazar is an attractive target due to the evidence of periodic oscillations observed at different wavelengths, with a solid proof of a 2.2-year modulation detected in the gamma -ray, UV, and optical bands. We aim to investigate the variability pattern of the X-ray emission using a more than 10-year-long light curve in order to robustly assess the presence or lack of a periodic behavior, evidence of which is only marginal so far. We conducted detailed statistical analyses, studying in particular the variability properties of the X-ray emission of by computing the Lomb-Scargle periodograms, which are suited for the analyses of unevenly sampled time series, and adopting epoch-folding techniques. We find a modulation pattern in the X-ray light curve of PG 1553+113 with a period of sim 1.4 years, which is about 35<!PCT!> shorter than the one observed in the gamma -ray domain. Our finding is in agreement with the recent spectro-polarimetric analyses and supports the presence of more dynamical phenomena simultaneously at work in the central engine of this quasar.

V0405 Dra: a new deep and low mass ratio contact binary with extremely fast decrease in the orbital period

Abstract V0405 Dra is a W UMa-type binary star. Basing on the TESS data, we have conducted an orbital period study and performed a light curve analysis for the system. The orbital period study reveals that the $O-C$ curve for V0405 Dra exhibits secular decrease at an extremely high rate of $dP/dt = -2.71 \times 10^{-6} day\cdot year^{-1}$, along with periodic variations characterized by an amplitude of $A_3 = 0.0032$ days and a period of $P_3 = 1.413 $ years. The orbital periodic change is possibly due to the light-travel time effect resulting from an additional third body in the system, for which we estimate a minimum mass of $M_3 = 0.77 M_{\odot}$. By employing the \add{2013 version of the} Wilson-Devinney (W-D \add{2013}) method to synthesize light curve, we derived photometric solutions indicating that V0405 Dra is a new deep ($f = 68.7\%$)and low-mass ratio ($q =0.175$) contact binary. The fast decrease in its orbital period is likely caused by mass transfer from the more massive primary star to the less massive secondary star, or due to angular momentum loss, which have been elaborated upon in the last section. With further mass transfer and loss of angular momentum, the binary will gradually evolve into a tighter contact configuration, eventually leading to a merge into a single star, following the evolutionary paths suggested for such deep and low mass ratio contact binaries.

Multicolour photometry of LEO mega-constellations Starlink and OneWeb

ABSTRACT The development of low earth orbit (LEO) mega-constellation fundamentally threatens ground-based optical astronomical observations. To study the photometric properties of the LEO mega-constellations, we used the Xinglong 50 cm telescope to conduct a large-sample, high-precision, and multicolour target-tracking photometry of two typical LEO mega-constellations: Starlink and OneWeb. Over a three-month observation period starting on 2022 January 1, we collected 1447 light curves of 404 satellites in four typical versions: Starlink v1.0, DarkSat, VisorSat, Starlink v1.5, and OneWeb. According to data statistics, Starlink v1.0 has the smallest median magnitude at clear and Sloan Digital Sky Survey gri band, and OneWeb is the dimmest bus. The brightness of Starlink v1.5 is slightly brighter than VisorSat. We construct a detailed photometric model with solar phase angle variations by calculating the illumination-visibility geometry based on the orbital parameters. Our data analysis shows that the solar phase angle is the significant characteristic which influencing Starlink satellites’ brightness, but it is not sensitive to OneWeb satellites. VisorSat and Starlink v1.5 versions, which are equipped with deployable visors, have significantly reduced scattered light compared to the previous Starlink v1.0 version. The multiband LOWESS and colour index are analysed in characterizing the energy and colour features of LEO mega-constellation satellites. This work found that the proportion of scattered sunlight mitigation achieved with VisorSat and Starlink v1.5 was 55.1 and 40.4 per cent, respectively. The colour index of different buses shows an evident clustering feature. Our observation and analysis could provide valuable quantitative data and photometric models, which can contribute to assessing the impact of LEO mega-constellations on astronomical observations.

Resolving the vicinity of supermassive black holes with gravitational microlensing

Abstract Upcoming wide field surveys will discover thousands of new strongly lensed quasars which will be monitored with unprecedented cadence by the Legacy Survey of Space and Time (LSST). Many of these quasars will undergo caustic-crossing events over the 10-year LSST survey, during which the quasar’s inner accretion disc crosses a caustic feature produced by an ensemble of microlenses. Such caustic-crossing events offer the unique opportunity to probe the vicinity of the central supermassive black hole, especially when combined with high cadence, multi-instrument follow-up triggered by LSST monitoring. To simulate the high cadence optical monitoring of caustic-crossing events, we use relativistic accretion disc models which leads to strong asymmetric features. We develop analysis methods to measure the Innermost Stable Circular Orbit (ISCO) crossing time of isolated caustic-crossing events and benchmark their performance on our simulations. We also use our simulations to train a Convolutional Neural Network (CNN) to infer the black hole mass, inclination angle, and impact angle directly from these light curves. As a pilot application of our methods, we used archival caustic-crossings of QSO 2237+0305 to estimate the black hole mass and inclination angle. From these data, two of these methods called the second derivative and wavelet methods measure an ISCO crossing time of 48.5 and 49.5 days, corresponding to a Kerr black hole mass of MBH = (1.5 ± 1.2) × 109M⊙ and MBH = (1.5 ± 1.3) × 109M⊙ respectively. The CNN inferred log10(MBH/M⊙) = 8.35 ± 0.30 when trained on Schwarzschild black hole simulations, and a moderate inclination of i = 45 ± 23○. These measurements are found to be consistent with previous estimates.

Toward an Auditory Virtual Observatory

The large ecosystem of observations generated by major space telescope missions can be remotely analyzed using interoperable virtual observatory technologies. In this context of astronomical big data analysis, sonification has the potential of adding a complementary dimension to visualization, enhancing the accessibility of the archives, and offering an alternative strategy to be used when overlapping issues are found in purely graphical representations. This article presents a collection of sonification and musification prototypes that explore the case studies of the MILES and STELIB stellar libraries from the Spanish Virtual Observatory and the Kepler and TESS light curve databases from the Space Telescope Science Institute archive. Using automation and deep learning algorithms, it offers a “palette” of resources that could be used in future developments oriented toward an auditory virtual observatory proposal. The work includes a user study with quantitative and qualitative feedback from specialized and nonspecialized users analyzing the use of sine waves and musical instrument mappings for revealing overlapped lines in galaxy transmission spectra, confirming the need for training and prior knowledge for the correct interpretation of accurate sonifications, and providing potential guidelines to inspire future designs of widely accepted auditory representations for outreach purposes.

An Online Framework for Fitting Fast Transient Lightcurves

Abstract The identification of extragalactic fast optical transients (eFOTs) as potential multi-messenger sources is one of the main challenges in time-domain astronomy. However, recent developments have allowed for probes of rapidly-evolving transients. With the increasing number of alert streams from optical time-domain surveys, the next paradigm is building technologies to rapidly identify the most interesting transients for follow-up. One effort to make this possible is the fitting of objects to a variety of eFOT lightcurve models such as kilonovae and γ-ray burst (GRB) afterglows. In this work, we describe a new framework designed to efficiently fit transients to light curve models and flag them for further follow-up. We describe the pipeline’s workflow and a handful of performance metrics, including the nominal sampling time for each model. We highlight as examples ZTF20abwysqy, the shortest long gamma ray burst discovered to date, and ZTF21abotose, a core-collapse supernova initially identified as a potential kilonova candidate.

What drives the variability in AGN? Explaining the UV-Xray disconnect through propagating fluctuations

ABSTRACT Intensive broad-band reverberation mapping campaigns have shown that AGN variability is significantly more complex than expected from disc reverberation of the variable X-ray illumination. The UV/optical variability is highly correlated and lagged, with longer lags at longer wavelengths as predicted, but the observed time-scales are longer than expected. Worse, the UV/optical light curves are not well correlated with the X-rays, which should drive them. Instead, we consider an intrinsically variable accretion disc, where slow mass accretion rate fluctuations are generated in the optical-UV disc, propagating down to modulate intrinsically faster X-ray variability from the central regions. We match our model to Fairall 9, a well-studied AGN with L ∼ 0.1LEdd, where the spectrum is dominated by the UV/EUV. Our model produces light curves where the X-rays and UV have very different fast variability, yet are well correlated on longer time-scales, as observed. It predicts that the intrinsic variability has optical/UV leading the X-rays, but including reverberation of the variable EUV from an inner wind produces a lagged bound-free continuum that matches the observed UV-optical lags. We conclude that optical/UV AGN variability is likely driven by intrinsic fluctuations within the disc, not X-ray reprocessing: the observed longer than expected lags are produced by reverberation of the EUV illuminating a wind, not by X-ray illumination of the disc: the increasing lag with increasing wavelength is produced by the increased contribution of the (constant lag) bound-free continuum to the spectrum, rather than indicating intrinsically larger reverberation distances for longer wavelengths.

Real-time Abnormal Detection of GWAC Light Curve based on Wavelet Transform Combined with GRU-Attention

Nowadays, astronomy has entered the era of Time-Domain Astronomy, and the study of the time-varying light curves of various types of objects is of great significance in revealing the physical properties and evolutionary history of celestial bodies. The Ground-based Wide Angle Cameras telescope, on which this paper is based, has observed more than 10 million light curves, and the detection of anomalies in the light curves can be used to rapidly detect transient rare phenomena such as microgravity lensing events from the massive data. However, the traditional statistically based anomaly detection methods cannot realize the fast processing of massive data. In this paper, we propose a Discrete Wavelet (DW)-Gate Recurrent Unit-Attention (GRU-Attention) light curve warning model. Wavelet transform has good effect on data noise reduction processing and feature extraction, which can provide richer and more stable input features for a neural network, and the neural network can provide more flexible and powerful output model for wavelet transform. Comparison experiments show an average improvement of 61% compared to the previous pure long-short-term memory unit (LSTM) model, and an average improvement of 53.5% compared to the previous GRU model. The efficiency and accuracy of anomaly detection in previous paper work are not good enough, the method proposed in this paper possesses higher efficiency and accuracy, which incorporates the Attention mechanism to find out the key parts of the light curve that determine the anomalies. These parts are assigned higher weights, and in the actual anomaly detection, the star is detected with 83.35% anomalies on average, and the DW-GRU-Attention model is compared with the DW-LSTM model, and the detection result f1 is improved by 5.75% on average, while having less training time, thus providing valuable information and guidance for astronomical observation and research.

A Comparative Study of TESS Light-curve Extraction Methods Applied to Blazars

Blazars are characterized by largely aperiodic variability on timescales ranging from minutes to decades across the electromagnetic spectrum. The Transiting Exoplanet Survey Satellite (TESS) mission provides continuous sampling of blazar variability on timescales ranging from tens of minutes to 27 days for a single sector observation. Proper removal of the background, thermal ramping, and onboard systematic effects are crucial to the extraction of a reliable blazar light curve. Multiple publicly available procedures have been created to correct for these effects. Using ground-based observations from the Zwicky Transit Facility (ZTF) and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as “ground truth” observations, we compare six different methods (regression, cotrending basis vectors, pixel-level decorrelation, eleanor, quaver, and simple differential photometry) to each other, and to our “ground truth” observations, to identify which methods properly correct light curves of a sample of 11 bright blazars, including the suspected neutrino source TXS 0506+056. In addition to comparing the resulting light curves, we compare the slopes of the power spectral densities, perform least-square fitting to simultaneous ZTF data, and quantify other statistical qualities generated from the light curves of each method. We find that only three of the six methods compared (simple differential photometry, eleanor, and quaver) produce TESS light curves consistent with the ground-based ZTF and ATLAS observations.

Kilonova Parameter Estimation with LSST at Vera C. Rubin Observatory

The upcoming Vera Rubin Observatory’s Legacy Survey of Space and Time (LSST) opens a new opportunity to rapidly survey the southern sky at optical wavelengths (i.e., ugrizy bands). In this study, we aim to test the possibility of using LSST observations to constrain the mass and velocity of different kilonova (KN) ejecta components from the observation of a combined set of light curves from afterglows of γ-ray bursts and KNe. We used a sample of simulated light curves from the aforementioned events as they would have been seen during the LSST survey to study how the choice of observing strategies impacts the parameter estimation. We found that the design of observing strategy that is the best compromise between light-curve coverage, observed filters, and reliability of the fit involves a high number of visits with long-gap pairs of about 4 hr every two nights in the same or different filters. The features of the observing strategy will allow us to recognize the different stages of the evolution of the light curve and gather observations in at least three filters.

Dramatic Rebrightening of the Type-changing Stripped-envelope Supernova SN 2023aew

Multipeaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility. Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (−17.4 mag) and long (∼100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew’s spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak’s photospheric phase, before Hα resurfaces again during the nebular phase. The nebular lines ([O i], [Ca ii], Mg i], Hα) exhibit a double-peaked structure that hints toward a clumpy or nonspherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than that of normal SESNe as well as requiring a very high 56Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak and also powers the second peak of the light curve through interaction of the SN with the circumstellar medium.

Discovery of a Second Eclipsing, ​​​​​Bursting Neutron Star Low-mass X-Ray Binary in the Globular Cluster Terzan 6

We have analyzed Chandra and Suzaku observations of the globular cluster Terzan 6, made when the recurrent transient GRS 1747–312 was in quiescence. Our analysis reveals the presence of a second eclipsing, bursting neutron star low-mass X-ray binary in the central regions of the cluster, in addition to GRS 1747–312. The new source, which we name Terzan 6 X2, is located only ∼0.″7 away from GRS 1747–312 in the 2021 Chandra images. The detection of a 5.14 ks long eclipse in the light curve of X2 at a time not predicted by the ephemeris of GRS 1747–312 confirms that it is an unrelated source. Using the Suzaku light curve from 2009, which, in addition to a type I X-ray burst, also showed an eclipse-like feature, we constrain the orbital period to be longer than 16.27 hr. The 0.5–10 keV luminosities of X2 vary in the range of ∼0.24–5.9 × 1034 erg s−1 on timescales of months to years. We have identified a plausible optical counterpart of X2 in Hubble Space Telescope F606W and F814W images. This star varied by 2.7 mag in V 606 between epochs separated by years. In the cluster color–magnitude diagram, the variable counterpart lies in the blue-straggler region when it was optically bright, about 1.1–1.7 mag above the main-sequence turn-off. From the orbital period–density relation of Roche-lobe filling stars, we find the mass-donor radius to be ≳0.8 R ⊙.

The Effects of r-Process Enrichment in Hydrogen-rich Supernovae

Core-collapse supernovae (SNe) are candidate sites for rapid neutron capture process (r-process) nucleosynthesis. We explore the effects of enrichment from r-process nuclei on the light curves of hydrogen-rich SNe and assess the detectability of these signatures. We modify the radiation hydrodynamics code, SuperNova Explosion Code, to include the approximate effects of opacity and radioactive heating from r-process elements in the supernova (SN) ejecta. We present models spanning a range of total r-process masses M r and their assumed radial distribution within the ejecta, finding that M r ≳ 10−2 M ⊙ is sufficient to induce appreciable differences in their light curves as compared to ordinary hydrogen-rich SNe (without any r-process elements). The primary photometric signatures of r-process enrichment include a shortening of the plateau phase, coinciding with the hydrogen-recombination photosphere retreating to the r-process-enriched layers, and a steeper post-plateau decline associated with a reddening of the SN colors. We compare our r-process-enriched models to ordinary SNe models and observational data, showing that yields of M r ≳ 10−2 M ⊙ are potentially detectable across several of the metrics used by transient observers, provided that r-process-rich layers are mixed at least halfway to the ejecta surface. This detectability threshold can roughly be reproduced analytically using a two-zone (kilonova-within-an-SN) picture. Assuming that a small fraction of SNe produce a detectable r-process yield of M r ≳ 10−2 M ⊙, and respecting constraints on the total Galactic production rate, we estimate that ≳103–104 SNe need be observed to find one r-enriched event, a feat that may become possible with the Vera Rubin Observatory.

3D hydrodynamic simulations of massive main-sequence stars II. Convective excitation and spectra of internal gravity waves

Abstract Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light curve observations. We present the oscillation properties of high resolution hydrodynamic simulations of a $25\;\mathrm{\mathrm{M}_\odot }$ star performed with the PPMstarcode. The model star has a convective core mass of $\approx 12\;\mathrm{\mathrm{M}_\odot }$ and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time scale in addition to significant power across a forest of low and high angular degree l modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes on the order of 10s of days. Thanks to the still significant power at higher l and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt–Väisälä frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary.